Dynamical scaling of imbibition in columnar geometries 
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Recent experiments of imbibition in columnar geometries show interfacial fluctuations whose 
dynamic scaling is not compatible with the usual non local model governed by surface tension 
that results from a macroscopic description. To explore this discrepancy, we exhaustively analyze 
numerical integrations of a phase-field model with dichotomic columnar disorder. We find that two 
distinct behaviors are possible depending on the capillary contrast between both values of disorder. 
In a high contrast case, where interface evolution is mainly dominated by the disorder, an inherent 
anomalous scaling is always observed. Moreover, in agreement with experimental work, the interface 
motion has to be described through a local model. On the other hand, in a lower contrast case, 
interface is dominated by interfacial tension and can be well modeled by a non local model. We 
have studied both spontaneous and forced-flow imbibition situations, giving a complete set of scaling 
exponent in each case, as well as, a comparison to the experimental results. 

PACS numbers: 47.56.-|-r, 68.35.Ct, 05.40.-a 



I. INTRODUCTION 

Fluid-fluid displacements in porous media is a sub- 
ject of much interest in industrial processes and material 
characterization so as in environmental problems rang- 
ing from petroleum recovery and irrigation to retention 
of waste waters [l|, 01 • We restrict our analysis to the case 
in which the invading fluid, that wets preferentially the 
medium, is more viscous than the displaced resident fluid 
giving rise to compact rough interfaces. This process 
can take place either spontaneously at constant pressure, 
spontaneous imbibition, or by the application of an ex- 
ternal pressure at constant injection rate, forced-flow im- 
bibition 2\. The surface can be represented by a single- 
valued function of position x and time t, h = h(x, t) 
Q. In many cases, when fluctuations have a thermal 
origin or can be reduced to that, the interfacial fluctu- 
ations are self-affine and follow the dynamic scaling of 
Family- Vicsek Then, a complete description of the 
dynamical growth is possible with only two independent 
exponents. However, when disorder is relevant enough 
to interfere with the geometry of the interface, a more 
generic scaling can apply, and one more independent ex- 
ponent is necessary to reach the dynamical description. 
It is the so-called anomalous scaling [5], and it has been 
observed in many different experimental and numerical 
situations during the last decade [1, 0, H, H, [13) [U • 

Several experiments of imbibition in distinct geome- 
tries have been proposed in the last years. There are 
experiments that use paper as a disordered medium, 
[II [H, [11, [H, [11] and others performed in Hele-Shaw 
cells (two parallel glass plates separated by a narrow dis- 
tance) with a random distribution of glass beads as a 
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disordered medium [l3, [iM [3| • Using the same geom- 
etry, other methods of generating a disordered medium 
have been explored, including random variations in gap 
spacing produced by a predesigned surface relief of the 
bottom plate [1, Hi [0, or by roughened plates [2l[. 
Focusing our attention on columnar geometries we have 
presented in previous works [l^ experimental stud- 
ies of forced-flow imbibition in a Hele-Shaw cell with 
a columnar quenched disorder, produced by dichotomic 
variations in the thickness of the bottom plate. We found 
that the interfacial dynamics followed an intrinsic anoma- 
lous scaling with varying exponents that were incompati- 
ble with the expected results from the usual macroscopic 
model, that gives characteristic exponents of a non lo- 
cal model. Instead, we obtained good agreement with an 
heuristic model of diffusively coupled columns present- 
ing local interactions. The essential change of behavior, 
which pass to be dominated by non local to local interac- 
tions, is obviously due to the persistence of the columnar 
disorder. However the detailed physical mechanism is not 
clear. 

In the present paper we use numerical integrations of a 
phase-field model with columnar disorder to explore this 
behavior. By imposing the columnar geometry of the in- 
terface motion at the macroscopic model, we analytically 
derive the heuristic local equation presented in Ref. 0. 
In addition, we show how both local and non local be- 
haviors are observed in the numerical model through the 
variation of a parameter modeling the contrast between 
the different capillary values present in the system. The 
outline of the paper is as follows. In Sec. [H] we intro- 
duce the phase-field model and the interfacial equations 
obtained in the sharp interface limit. In Sec. Illll we show 
the notions of generic scaling used to characterize the in- 
terfacial dynamics. Section HVl is devoted to analyze the 
numerical results in the low capillary contrast case for 
both spontaneous and forced-flow imbibition situations. 
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Section |V] deals with the case of high capillary contrast 
and its connection with local growth models. Finally, in 
Sec.|Vl]we discuss the physical contents of phase- field and 
macroscopic models as well as their relevance to explain 
the experimental results. 



II. THE PHASE-FIELD AND MACROSCOPIC 
MODEL 

Our numerical results shall be performed using the so- 
called phase-field model [13, [H, HJ]. Such a model is 
based on the introduction of an order parameter which 
can take two limit values, = ±(/)e representing the 
two phases liquid/air of the system. The phase-field dy- 
namics is controlled by a conserved equation based on a 
Ginzburg-Landau formulation, dcjj/dt — VMV^, where 
fi = SJ- /S(f> is the chemical potential and the free energy 
IS given by the functional T[(l>] = J dr[V{(l>) + i(eV0)2]. 
The phase-field equation then reads as 



dt 



VA/V 



(1) 



where M is a parameter that is taken constant in the 
liquid phase and zero in the air phase. V{(j)) is a poten- 
tial taken as V{4>) — — -f — rj{r)(j) and defines 
two stable phases through the double well potential; the 
destabilizing linear term accounts the effect of a capil- 
lary force that makes the interface to advance. The ef- 
fect of an inhomogeneous capillarity is added by using a 
dichotomic capillary noise with the values: 



no 



(2) 



We consider a columnar disorder r]{x) defined by single 
tracks of lateral size Ld distributed along the x-direction 
in such a way that tracks with the high disorder value 
V = Vo/i^ — Va) occupy 35% of the system length L. 
This is the same kind of disorder reported in the exper- 
imental work of Ref. 0]. Note that, for a given disorder 
realization, tracks wider than Ld are obtained when two 
or more unit tracks are placed adjacently. 



A. The macroscopic description of imbibition 

The use of phase-field models to reproduce imbibition 
experiments is based on the ability to get the same results 
than those obtained from a macroscopic model. Indeed, 
in the so-called sharp interface limit e — > 0, a matched 
asymptotic expansion of the field (j) around a kink solu- 
tion of Eq. ([1]), 00 = —4'eqtcLnh{w/ ^/{2)), can be per- 
formed, recovering then the basic macroscopic equations 
for the usual pressure p, velocity v, interfacial curvature 
K and the columnar capillary disorder which we called 



V^p = 0, 

Apint ^ <Tfi + r]{x). 



(3) 
(4) 
(5) 



The macroscopic variables and parameters are defined 
from the phase- field formulation as [2^ : 



K 



M 
Wo 



dw 



\dw) ' 



fii being the first order term of the expansion on e 
of the chemical potential and a playing the role of a 
interfacial tension. The three equations of the macro- 
scopic model are well known from phenomenological 
arguments involving conservation laws. Darcy's law, Eq. 

arises from an averaging procedure of Navier-Stokes 
equations at low Reynolds number, when the geometry 
of the Hclc-shaw cell is imposed. The Laplace equation, 
Eq. (|4]), comes from imposing incompressibility of the 
liquid, and the Gibbs-Thomson relation, Eq. ([5]), comes 
from a principle of minimum interfacial energy. The 
capillary pressure at the interface can be expressed as 
ri{x) ~ 2acos9/b, where 6 corresponds to the contact 
angle, and b is the distance between the plates of 
Hele-Shaw cell. 

In the experimental work reported in Ref. 9], the ran- 
dom gap distribution is constructed by using a fiberglass 
substrate with a pattern of copper tracks attached to 
the bottom plate of a Hele-Shaw cell. Tracks have a 
thickness d with a lateral size of Ld and are distributed 
along the lateral direction x without overlap. Therefore, 
the gap of the Hele-Shaw cell has a dichotomic variation 
with two possible values b and b — d. Although in that 
case the capillarity is a 3-D effect of the cell, we can 
relate our numerical parameters of Eq. ^ as Tyo = 1/& 
and TjA = d/b. Note that the parameter rjA is related to 
the capillary contrast between both values of disorder. 



B. Spontaneous and forced-flow imbibition 

In our study, we will consider both situations of spon- 
taneous and forced-flow imbibition by choosing con- 
veniently the boundary conditions into the phase-field 
model [26l |. For spontaneous imbibition an applied con- 
stant pressure is imposed at the origin of the cell fi{x, y = 
0) = ^La- In contrast, a pressure gradient has to be im- 
posed at the origin, Kdy^\y=Q = —Vm, to reproduce 
forced-flow imbibition. Therefore, the main difference 
between both cases is found on the mean velocity of 
the interface. While in the forced-flow case, the inter- 
face evolves with the imposed constant velocity Kn, in 
the spontaneous imbibition case, the averaged interfa- 
cial height H{t) follows the so-called Washburn's law 
H{t) ~ t^/"^. An exact expression of such an evolution 



3 



can be obtained solving the following equation: 
dH Kifia + iv)) 



tension term aKIkr and the drift term \k\Vm \22 



dt 



Hit) 



(6) 



that comes from ([3]), (|3|), and (O with k = Q 
and (rj) = {iri{x))x- The expression for H reads 
H{t) = ^i7(0)2 + 2at where a = K{{tj) + fia)- 
Therefore, spontaneous imbibition has a slowing-down 
dynamics with a mean velocity Vm{t) ^ t~^^'^. 




(9) 



In the forced-flow case, Vm is constant and thus the 
crossover length is just a static length scale separating 
two different regimes. However, as it has been observed in 
several numerical results [H,|2a,[13| , this crossover length 
acts as a cutoff for the interface fluctuations growth due 
to the interface being asymptotically flat on length scales 
larger than ^ x ■ 



C. Equations for the interface 



From the macroscopic equations, Eqns. ([3]), ([4]) and 
([5]), it is possible to obtain an equation for the moving 
interface by means of a Green analysis. The Green func- 
tion G{r, r') in our case follows a Poisson equation with a 
unit source at r', V'^G{r/r') = 5{r — r') evaluated at the 
plane 2-D. The expression for the interface is obtained 
using the Green identity [131 : 

dr'[p{r')V'^G{r/r') - G{r /r')V'^p{r')] = 
J^^ ds' ■ pis')V'G{s, s') - ds' ■ Gis, s')V'p{s'), 

(7) 

integrated over the volume of the liquid Ql = {x,0 < 
y < h{x,t)}, h{x,t) being the interfacial position. De- 
pending on the intensity of the quenched noise two cases 
are relevant in our study. 



1. Low capillary contrast. Linearized equations 

We assume that the quenched noise does not impose 
any special geometry and the interface can be linearized 
around its averaged value, h{x, t) = H{t) + Sh{x, t), with 
H{t) — {h{x,t))j.. Keeping the first order on 5h fluc- 
tuations in Eq. ([7]), and imposing Darcy's law, Eq. ([3]), 
at the interface boundary conditioii, one can obtain the 
linearized interface equation 0, [23 , [26j , which in Fourier 
space reads as 

5hk = -(jK\k\k^5hk - Vm\k\5hk + K\k\r]k, (8) 

where we have supposed that correlations do not grow 
up faster in time than the mean height of the inter- 
face, \k\H{t) ^ 1. Under this limit, interface fluctua- 
tions follow the same Eg. (l8l) both for spontaneous and 
forced- flow imbibition [2J, [lo, HI] ■ However, it is worth 
to mention here that the presence of Washburn's law 
{Vm{t) t^^/^) in spontaneous imbibition give rise to 
dynamic crossover lengths and therefore, as it has been 
pointed out in Ref. [29|, a rich variety of different scal- 
ing regimes can be observed. The crossover length scale 
can be explicitly seen as a balance between the surface 
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FIG. 1: Scheme of the interface hi{t) advancing through an 
effective channel Li composed of several unit tracks Ld- The 
gray stripes correspond to tracks with a high capillary noise, 
V = Vo/{^ — Va), and the white ones to tracks with a low 
capillary noise, rj = rjo- 



2. High capillary contrast. Coupled channel equations 

The situation becomes quite different when the capil- 
lary contrast is increased. In the experimental work of 
Ref. [§], where forced-flow imbibition is studied for high 
capillary contrast, it is observed that the interface mo- 
tion can be modeled through a phenomenological local 
equation. Our purpose now is to derive such an equa- 
tion directly from the macroscopic model. To do this, we 
shall assume that the noise is so strong that the interface 
adopts the columnar geometry of the disorder. The pro- 
cedure consists then on integrating Eq. ([7]) in a closed sur- 
face along the proflle of the i-th effective channel. Here, 
an effective channel, or simply channel i, is defined as a 
group of single adjacent tracks of which the majority has 
the same disorder value, in such a way that the interface 
advances as a compact surface through each channel. An 
example is depicted in Fig. [1] where the channel is com- 
posed of several tracks of which the majority has a high 
disorder value (gray tracks) and therefore, the averaged 
disorder of the channel rji is larger than the mean disor- 
der of the whole system (77). Note that the surrounding 
channels i — 1 and i + 1 must have a mean disorder value 
rji < (r/). The width of the channel is defined as L^. A 
numerical example of these channels can be seen in Fig. 
2. 




(a)Forced-flow imbibition (b)Spontancous imbibition 



FIG. 2: Interface profiles at equal times but different values of capillary contrast rjA- The columnar disorder is also plotted: 
gray tracks are points where the dichotomic disorder takes its highest value of Eq. 



We are considering the forced-flow imbibition case with 
the following boundary conditions at the top and the 
bottom of the channel: 



dyPix,y)\y^a ^ 

K 



-iv), 



(10) 

(11) 

(12) 



where the pressure at the origin has been estimated by 
using the relation = — ^ , which comes essen- 

tially from Eq. pd|) , taking that Pint = — {r]) is the mean 
capillary pressure of the whole system. It means that 
the pressure at the origin is changing on time in order 
to get a mean constant velocity for the whole interface. 
Using these boundary conditions, and taking the general 
expression for the two-dimensional Green function 



G{x -x',y- y') 



dfc- 



h.2 I h.2 



we can evaluate the different boundary integrals of the 
top and bottom segments of Eq. ([7]) as 



Xi-l 
Xi-1 



dx n ■ [p{x,y)VG{x,y)\y^j^^ 

1 

da: n • [G{x,y)Vp{x,y)\y^j^^ 

1 

dx n ■ [p{x,y)VG{x,y)\y^^ 

1 

da; n- [G{x,y)Vp{x,y)\y^^^ 







ttK ^ 




VrnHit) 


iv) 


2K 


2 ' 






2K 


nK 



where C — J^du^^^^^^^, with X = Li/L being a cut off 
due to the finite size L of the system. We have also sup- 
posed enough wide channels to ensure that Li > hi(t), 
which means that we are taking the initial times, before 



interface gets saturated. The nomenclature used for the 
Green function means that it is evaluated at the inter- 
face, G{x,y) = G{0,hi/x,y). Therefore, Eq. ^ can be 
written as 

-p(0, h,) = j^dy' \p{x, y')dxG{x, y') - G{x, y')dxp{x, y) 

hi 

<iy' P{x, y')dxG{x, y') - G{x, y')dxp{x, y) 



iv) 

-a.ihi{t) + aiVm —■ 



(13) 



The two first terms of the right-hand side are due to the 
flow between neighbor channels. We have defined the 
parameter Ui = In order to get an equation for 

the time evolution of the interface hi{t), we define the 
following coupling coefficients between channels as the 
ratio between channel flow and height differences 



A =- 



J°dy' p{x, y')d^G{x, y') - G{x, y')d^p{x, y'] 



A-i 



ai{ht+i - hi) 
Jl^dy' _p(x, y')d,^G{x, y') - G{x, y')d^p{x, y'] 



ai{hi 



The coupling variable Di has to be understood as a dif- 
fusion coefficient which depends on each channel i. As 
a general case, it may be taken as a random variable. 
Moreover, we shall assume that Di does not vary in time 
during the initial times, before interface gets saturated. 
Then, we can write Eq. (|13p as an inhomogeneous diffu- 
sion equation between channels: 



-p.(t)-(77) . (14) 



hi{t) = VAV/li +Vrr, + ^ 

iQi 

where Pi{t) = p(0, hi) corresponds to the pressure at the 
interface. We are assuming that the pressure at the in- 
terface is time-dependent, which is based on the experi- 
mental results reported in Ref . [oj . In such experimental 
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work, the local velocity of the interface at each chan- 
nel follows an expression similar to Washburn's law un- 
til it reaches the saturation value Vm- In order to take 
into account such behavior into the equation, we con- 
sider that the pressure at the interface can be expressed 
as Pi{t) ^ —Apint — cui(t), where the pressure difference 
Apint is given by the usual Gibbs- Thomson relation, Eq. 
(5), taking a negligible atmospheric pressure. The term 
cui{t) is an effective kinetic term due to the local capillary 
forces at each channel, being c an arbritrary constant, 
and can be explained in terms of the mass conservation. 
When the interface goes through a channel j of a high 
capillary disorder {rji > {-q}), its local velocity tends to 
initially increase up to a nominal value. In contrast, since 
we are imposing a constant velocity for the whole inter- 
face, the local velocity at the neighbor channel i + 1 with 
a lower capillary disorder tends to decrease down to a 
nominal value. After reaching such nominal value in both 
cases, the local velocity decrease or increase asymptoti- 
cally to the saturation value following the Washburn's be- 
havior due to the capillary forces of each channel. There- 
fore, we are taking that Ui{t) i^^i^^^^ for t > 0, where 
is a random variable defined as Ci = {Vi ~ (^))/l'7i ~ ('7)li 
which takes the values £,i = +1 at the channel with the 
highest capillary value (r]i > (rj)), and = — 1 at the 
channel with the lowest capillary value (ry^ < (77)). In 
addition, we also suppose that the curvature of the in- 
terface can be approximated as a constant value kq, and 
only its sign depends on each channel as Ki ~ ^i^o- 
Then, rewriting the last term on the right-hand side of 
Eq. (HH), we get the final expression 

h,{t) = VD.Vh, + V.,n+^^ {v, + a,r 1/2) , (15) 

where we have defined the new constants as = {\'ni ^ 
(77)1 -I- i7/vo)/2ai, and = CFc/2ai. 

We thus conclude that when the capillary disorder is large 
enough, the columnar geometry of the system leads to a 
local description for the interface motion. 



which it is related to the others exponents through the 
scaling relation /? = a/ z. To study local growth, we shall 
define local widths averaged on windows of size £ < L, 
w{i,t) — {{h{x,t) — h{x,t)^y^y/'^. Local growth is also 
given as a power law w{£,t) ~ ic{t)°' for £ > £c and 
w{£,t) - 4(i)"-"'- for £ < 4, where aioc is the local 
roughness exponent. In the case of self-afline growth, 
global scaling coincides with local scaling, aioc = ct, and 
the fluctuating interface is well characterized by only two 
independent exponents. However, following the method 
shown in Ref. [5], we must use at least three indepen- 
dent exponents in a more general case. In our study we 
use time variation of the local width for several window 
lengths: w{£,t) ~ g{£t-^/^), with 

{u"'°<=,u < 1 
(16) 
const, u ^ 1 

obtaining direct measures of the exponents and (3* = 
3 — s^ioc^ Indirect measures of z and a can be ob- 

" z 

tained through a collapse of the individual figures as 
w{£,t)£~" ^ g'{£t~^/^), where now the scaling function 
goes as for u << 1 and u~" for u >> 1. Fi- 

nally, although the collapse gives an indirect measure of 
exponents in a very robust way it is worth to have a di- 
rect measure of at least three exponents. Thus we also 
use the evolution of the power spectral density S{k, t) = 
(|ft,(fc,t)p) that scales as S{k,t) = j^hrcs{kt^^^), with 

( + < 1 

sin) = <^ (17) 

Depicting S{k, t) at different times we have the direct 
measure of the spectral exponent as that coincides ei- 
ther with the global roughness exponent a, when power 
spectrum do not shift in time (self-afHne or superrough 
scalings) or with the local one a/oc, when a temporal shift 
is observed (intrinsic anomalous scaling). 



III. GENERIC SCALING LAWS AND THE 
CHARACTERIZATION OF THE INTERFACIAL 
DYNAMICS 

Rough interfaces grow exhibiting power laws in both 
a horizontal correlation length 4 ~ i^^^, that accounts 
the range of correlation, and a vertical growth length, 
like the interface width W{t,L) ~ . W is defined 
as the deviation of the height h{x,t) as W{L,t) = 
{{h{x,t) — h{x,t))^y^'^ (where (...) and — mean sample 
and spatial average respectively), z and /3 are the so- 
called dynamic and growth exponents, which completely 
characterize the growth of self-affine processes. The sat- 
uration of the surface occurs in a saturation time ts, 
when the correlation length reaches the system length 
£c{ts) — L. Above this time, interfacial width scales as 
W{t > ts,L) ~ L", a being the roughness exponent 



IV. LOW CAPILLARY CONTRAST 

In this section we deal with the case of low capillary 
contrast in both forced-flow and spontaneous imbibition 
cases. It is worth to mention here that the numerical 
parameters used in all numerical results have been e — 1, 
M — 1 (dimensionless units), and 770 = 0.3. Equation 
^ has been integrated over a system of lateral size L = 
256 using a spatial grid of Ax = 1.0 and a time step of 
Ai — 0.01. The minimum length of the track disorder 
has been taken as = 2. 

A. Forced-flow imbibition 

We start to study a regime of low capillary pressures. 
The used value for the capillary contrast has been of 
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FIG. 3: Statistical analysis of the interface fluctuations in the 
forced- flow case with a low capillary contrast, 77A = 0.66. (a) 
Local width w{£,t) evaluated at different window sizes. The 
global and local growth exponents can be measured directly 
from the data. The inset shows the best data collapse for 
the scaling function using the values of 2; = 3.0 and a — 1.5. 
It also suggests a local roughness exponent of a;oc ~ 1-0. 
(b) Interface power spectrum calculated at different times. It 
shows a roughness exponent of as = 1.3 ± 0.2. 



rjA — 0.66 and the mean velocity has been fixed to 
Vra — 0.0025. The shape of the interface for a given real- 
ization is shown in Fig. [^Ji., case -qA = 0.66. We can see a 
smooth interface vifhich is slightly correlated with the dis- 
order. The results of the roughness analysis are shown in 
Fig. [31 From the local width u'(^, t) computed at different 
window sizes we can obtain the global and local growth 
exponents, /3 = 0.49±0.05 and /3* = 0.13±0.05. The best 
collapse of theses curves (shown at the inset of Fig. [5^ is 
obtained tuning the values oi z = 3.0 and a. = 1.5. The 
slopes of the scaling function agree with the previously 
calculated exponents suggesting also a value for the local 
roughness exponent of aioc — 1 which is corroborated at 
the power spectrum calculated at different times. Since 
there is not any temporal shift between the lines of the 
power spectrum, we can assume that the interface fluc- 
tuations are described within the superrough anomalous 
scaling, and therefore aioc = 1. In addition, the spectral 
roughness exponent obtained from the power spectrum 
as = 1.3 ± 0.2 corresponds to the global roughness ex- 
ponent, which is in agreement with the value obtained 
previously. These measured exponents are compatible 



with those obtained by the linear Eq. ([8]) with a con- 
stant velocity Vm- Rescaling such a linear equation by 
the transformation x —> bx,t — > b^t,h — * b'^h, we have 
trivially that z = 3 and a = 1.5, and assuming super- 
roughness {aioc — l)j we get the remainder exponents 
a.s P — a/z — 0.5, P* — (3 — = 0.17, in agreement 
with the measured exponents. Therefore, we can con- 
clude that this regime is well modelled by the non local 
and linear equation ([8]), taking into account only the sur- 
face tension regime. 

If we increase the mean velocity of the interface, then the 
static crossover length, Eq. (9), decreases, obtaining that 
interface fluctuations saturate earlier, at the time when 
the correlation length £c reaches the crossover length. 



B. Spontaneous imbibition 

In spontaneous imbibition, the crossover length scale, 
Eq. ([9]), becomes a dynamical scale and different regimes 
can be observed depending on the velocity of the inter- 
face [lO]- For low velocities, the initial correlation length 
tc ~ t^^^ is below the crossover length meaning that the 
relevant mechanism to damping the interface fluctuations 
is the surface tension with the characteristic dynamical 
exponent of z = 3. On the other hand, for higher veloci- 
ties, the crossover length acts as an effective correlation 
length of the interface fluctuations, giving rise to the gen- 
uine exponent of z = 4. In order to study both regimes, 
we have controlled the initial velocity of the interface by 
choosing the initial height of the interface conveniently. 
We impose an initial height of H{Q) = 199 to study a low 
velocity regime, whereas a higher velocity regime will be 
achieved by simply putting H{Q) — 1. 



1. Low velocity regime 

Typical shapes of interfaces in spontaneous imbibition 
are depicted in Fig. ^Bp. In the low capillary regime, case 
of r]A = 0.64, the interface is weakly correlated with dis- 
order. A roughness analysis shown in Fig. 2] gives the 
exponents /3 = 0.5 ± 0.04 and (i* = 0.18 ± 0.04 from a 
direct measure of the growing local width. The best data 
collapse on these figures provides the exponents z = 3.0 
and a = 1.5. Likewise, the slope of the scaling func- 
tion gives aioc = 1, which is corroborated by the power 
spectrum shown in Fig. , where we get a spectral 
roughness exponent of a = 1.35 ± 0.2, without temporal 
shift between the curves, indicating we are dealing with a 
superrough scaling. Hence, the measured exponents are 
the same than those obtained in the forced-flow case. It 
was actually expected, since the relevant terms of Eq. ^ 
at low velocities are the same in both cases. 
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FIG. 4: Statistical analysis of the interface fluctuations in 
the spontaneous case within a low velocity regime, and with 
a low capillary contrast, rjA = 0.64. (a) Local width w{£, t) 
evaluated at different window sizes. The global and local 
growth exponents can be measured directly from the data. 
The inset shows the best data collapse for the scaling function 
using the values of 2 = 3.0 and a = 1.5. It also suggests a 
local roughness exponent of qjoc = 1.0. (b) Interface power 
spectrum calculated at different times. It shows a roughness 
exponent of Qs = 1.35 ± 0.2. 
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FIG. 5: Statistical analysis of the interface fluctuations in 
the spontaneous case within a high velocity regime, and with 
a low capillary contrast, rjA ~ 0.64. (a) Local width w{£, t) 
evaluated at different window sizes. The global and local 
growth exponents can be measured directly from the data. 
The inset shows the best data collapse for the scaling function 
using the values of 2 = 4.0 and a — 1.5. It also suggests a 
local roughness exponent of aioc = 1.0. (b) Interface power 
spectrum calculated at different times. It shows a roughness 
exponent of as = 1.35 ± 0.2. 



2. High velocity regime 

When the initial interface velocity is increased, the 
velocity-dependent term of Eq. ([5]) starts to be relevant 
and the new regime adopts the dynamics of the crossover 
length t^^^, getting then the expected dynami- 

cal exponent z — 4 f25|, keeping the spatial structure 
with the same roughness exponent as before, a = 1.5. 
The numerical results are presented in Fig. [5l obtaining 
/? = 0.37 ± 0.03, (3* = 0.12 ± 0.03 and a, = 1.35 ± 0.2 
from direct measures, and a = 1.5, z = 4.0, a/oc — 1 
from the data collapse of the local widths, which are also 
in agreement with the linear description of Eq. ([8]). 



V. HIGH CAPILLARY CONTRAST 

When the capillary contrast is increased the shapes of 
interfaces become sharper in both cases of imbibition, as 
it can be seen in Fig. [21 The effect of disorder is so strong 



that a kind of columnar geometry is also imposed on the 
interface. As we will see, two important points of the dy- 
namics of interface fluctuations arise as effect of increas- 
ing the capillary contrast. First, interface motion seems 
to be described by local effects instead of the non local 
effects observed in the case of low capillary contrast. Sec- 
ond, scaling type changes from superrough to anomalous 
intrinsic. As in the case of low contrast, forced-flow and 
spontaneous imbibition cases present distinct patterns of 
fluctuations and they need to be studied separately. 

A. Forced-flow imbibition 

We start to study the case of a capillary contrast 
riA — 0.72 (see Fig. [2^). In Fig. there are plotted the 
local interface velocity at two different points x of the 
system (solid lines) and the mean velocity Vm (dashed 
line). The curve above the mean velocity corresponds to 
a channel i with a high disorder value, that is 77^ > (77), 
being rji the disorder of the channel i and (ry) the mean 
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FIG. 6: Interface velocity v{t) evaluated at two different 
points X of the system (solid lines) and the mean interface 
velocity Vm (dashed line) in both imbibition cases: (a) Forced- 
flow imbibition with Vm = 0.0025, and (b) spontaneous im- 
bibition within a high velocity regime. The curve above the 
mean velocity correspond to a channel with a high capillary 
disorder value (gray tracks in Fig. 2) and the curve below the 
mean velocity correspond to a channel with a low capillary 
disorder value (white tracks in Fig. 2). 




FIG. 7: Statistical analysis of the interface fluctuations in 
the forced- flow case with a high capillary contrast, rjA ~ 0.72. 
(a) Local width w{£, i) evaluated at different window sizes. 
The global and local growth exponents can be measured 
directly from the data. The inset shows the best data 
collapse for the scaling function using the values of z = 2.2 
and Q = 1.1. It also suggests a local roughness exponent of 
otioc = 0.7. (b) Interface power spectrum calculated at dif- 
ferent times. It shows a roughness exponent of = 0.65±0.2. 



disorder of the whole system. The curve below the mean 
velocity corresponds to the next channel i + 1 with a 
lower disorder value, that is 77^+1 < (77). We can see 
that both profiles can be locally described by Eq. 
v± ^ Vmi{v + dt~^/'^), where -I- and — mean the channel 
with rji — (77) > and rji — (rj) < respectively. Since 
the averaged velocity of the interface is low enough, the 
velocity w_ can be initially negative . 
As in the previous section we are interested in calculat- 
ing the different scaling exponents. The scaling analysis 
is shown in Fig. [71 From direct measure of the inter- 
facial local width slopes we get f3 = 0.52 ± 0.05, and 
/?* = 0.26 ± 0.06. The best data collapse is obtained us- 
ing a = 1.1 and z = 2.2. The power spectrum evolution 
changes with respect to the low contrast case since now 
there is a temporal shift between the curves, indicating 
the presence of inherent anomalous scaling. Hence, the 
spectral roughness exponent — 0.65 ± 0.2 must be in- 
terpreted as the local roughness exponent, which is in 
agreement with the slopes of the collapsed scaling func- 
tion. Note that these exponents are very close to the ex- 
perimental exponents reported in Ref. [01, /3 = 0.5±0.04, 
P* = 0.25 ± 0.03, a = 1.0 ± 0.1. As shown in this ref- 



erence, the shape of the interfaces and the analysis of 
fluctuations are well reproduced by Eq. (fTS]) . which takes 
into account strong diffusive coupled channels. There- 
fore, we can conclude that simple numerical integrations 
of the phase-field model reproduce both the shape and 
the scaling analysis of interfaces. 



1. From low to high capillary contrast 

When we increase even more the values of the capil- 
lary contrast, interfaces become more correlated to the 
columnar disorder presenting quantitative changes on the 
scaling behavior. For instance, in Fig. [8^ there is plotted 
the local width computed in a small window of length 
i = L/128. It allows us to calculate the local growth 
exponent (3* for different capillary contrasts. The com- 
plete set of scaling exponents is presented in Table ([l|. 
We can see that for high capillary contrasts {?]a > 0.7), 
interface fluctuations are always described by intrinsic 
anomalous scaling. In addition, the dynamics of the cor- 
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(a) 




FIG. 8: Interface local width calculated in a small window 
of size I — I//128 for different capillary contrasts. The local 
growth exponent /3* can be obtained directly from the data, 
(a) Forced-flow imbibition, Vm ~ 0.0025. (b) Spontaneous 
imbibition within a high velocity regime. 



length and Di ^ D can be taken as constant, obtaining 
then the regime described by z ~ 2. 

B. Spontaneous Imbibition 

In Fig. there is plotted four interface profiles eval- 
uated at the same time but different capillary contrast. 
The velocity profiles of the interface are shown in Fig. 
|BJd. Since now there is not any imposed velocity, the ve- 
locity of the interface follows the Washburn's law in each 
channel of different noise value. 

When we increase the parameter rjA a transition to a de- 
coupled state {(} = (}*= 0.5, z = oo) also appears in 
the case of spontaneous imbibition. However, there is an 
important difference with the forced-flow case. As be- 
fore, we calculate the local width in a small window size 
i = L/128 for different capillary contrasts (Fig.[5)3). We 
observe that now, the local growth exponent (3* changes 
suddenly to (3* = 0.5, indicating that the interface ad- 
vances completely decoupled. The transition to the de- 
coupled state can be now discontinuous in time for each 
value of the capillary contrast. It seems that there ex- 
ists a length (diVA) above which the fluctuations become 
decoupled. Therefore, for a given value of the capillary 
contrast, the clusters of interface with a size id will be- 
come decoupled between each other at the time td ~ 
when correlation length £c ^ t^^^ reaches the length id- 
Above td, the local description is not valid anymore. 



relation is subdiffusive, ranging from z = 2 to z ^ cxd in 
the highest contrast value. We must interprete this ex- 
treme case as having fluctuations completely decoupled. 
It is characterized by the same local and global growth 
exponent /3 = (3* = 0.5, which implies a dynamic expo- 
nent of z = oo, and it can be understood as a regime 
where the correlation length ic ~ t^^^ does not grow in 
time anymore. Experimentally, it has also been observed 
in the work carried out in Ref. [lol |. It is worth noting 
that there are two exponents /3 = 0.5 and aioc = 0.5 that 
remain constant. 

Alternatively to the phase-field model results, this phe- 
nomenon can be reproduced by Eg. (|15p taking a ran- 
dom diffusion coefficient D{x) [3l|. By choosing a prob- 
ability density = NaD-'' f,{D / Dmax), fc being a 
cut off function, Eq. (jl5p predicts a growth exponent 
(3 = aioc — 0.5 independently of the a value, a rough- 
ness exponent \ < a < oo, and a dynamic exponent 
2 < z < oo, which is in accordance with the numerical 
values obtained using the phase-field model (see Table 
m . In this sense, we can say that the diffusion coefficient 
Di of Eq. depends on the capillary contrast of the 
system. For very high capillary contrasts [r]A > 0.8) the 
coupling coefficient can be taken as ~ 0, obtaining 
then the decoupled state observed numerically. On the 
other hand, for lower capillary contrasts {rjA — 0.72) the 
variations of Di occur at scales larger than the correlation 



VI. CONCLUSIONS 

By means of numerical integrations of a phase-field 
model we find that there are strong differences between 
the dynamics of fluctuations in the cases of low and 
high capillary contrast with columnar disorder. Although 
these differences also exist in a case of quenched noise 
composed of squares |30|] , they are more dramatic when 
the quenched noise is of columnar type. It turns out that 
the persistence of the noise forces the interface to adopt 
the same geometry, changing the nature of the interface 
motion. 

In a low capillary contrast case, interfaces are superrough 
with a dynamics dominated by surface tension with ex- 
ponents z = 3 when the velocity is nearly constant or 
z — A when the velocity varies with Washburn's law. Fur- 
thermore, as interfaces are smooth and can be linearized 
around their mean value, a simple non local model for 
the interfacial evolution can be used to explain the ob- 
served dynamical scaling of fluctuations. On the other 
hand, when the capillary contrast is increased, interfaces 
are sharper and the correlation with disorder is more ev- 
ident. The observed dynamical scaling corresponds then 
to an anomalous scaling description with a clear tempo- 
ral shift at the power spectrum, and a subdiffusive be- 
havior with dynamical exponents ranging from z = 2 to 
z — > oo, depending on the strength of capillary forces. 
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TABLE I: Complete set of scaling exponents for forced-flow imbibition when the capillary contrast r)A is increased. 





a 




a. 


z 


a 


13* 


Scaling class 




0.66 


1.5 


1 


1.2 


3 


0.49 


0.12 


Superroughness 


Non-local description 


0.72 


1.1 


0.6 


0.6 


2.2 


0.52 


0.26 


Intrinsic A. 




0.74 


1.85 


0.55 


0.55 


3.7 


0.51 


0.35 


Intrinsic A. 


Local 


0.78 


2.5 


0.5 


0.5 


5.1 


0.5 


0.41 


Intrinsic A. 


Description 


0..S2 




0.5 


0.5 


3C', 


0.5 


0.5 


Docoupk'd state 





A prominent point to remark is that this behavior can 

be explained by a local model made of coupled channels 
with a fluctuating force following Darcy's law. One can 
interpret that in the high contrast case, the columnar 
disorder induces the existence of channels with more or 
less coupling, eliminating completely the non local char- 
acter of imbibition in homogeneous geometries. Finally, 
the difference between forced-flow and spontaneous im- 
bibition has also been elucidated in the high capillary 
contrast regime. While in the forced-flow case the inter- 
face gets completely decoupled above a critical capillary 



contrast, in spontaneous imbibition the same decoupled 
state can be suddenly achieved for a given value of cap- 
illary contrast. 
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